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GENERALIZED  CONICAL-FLOW  FIELDS  IN 
SUPERSONIC  WING  THEORY 
By  Harvard  Lomax  and  Max.  A.  Heaslet 


SUMMARY 


Linearized,  compressible— flow  analysis  is  applied  to  the  study  of 
quasi— conical  supersonic  wing  theory.  Single— integral  equations  are 
derived  which  relate  either  the  loading  to  the  shape  of  a  lifting 
surface  or  the  thickness  of  a  symmetrical  wing  to  the  pressure  distribu¬ 
tion  for  triangular  wings  with  subsonic  leading  edg9S.  The  forms  of 
these  equations  and  their  inversions  are  simplified  through  the  intro¬ 
duction  of  the  finite  part  and  the  generalized  principal  part  of  an 
integral. 

Applications  of  the  theory,  in  the  lifting  case,  include  previously 
known  results.  In  the  nonlifting  case,  it  is  shown  that  for  a  specified 
pressure  distribution  the  theory  does  not  always  predict  a  unique  thick¬ 
ness  distribution.  This  is  demonstrated  for  a  triangiilar  plan  form  hav¬ 
ing  a  constant  pressure  gradient  in  the  stream  direction. 

INTRODUCTION 


If  a  sufficiently  thin  wing  at  a  small  angle  of  attack  is  placed  in 
a  uniform  stream,  its  aerodynamic  properties  can  be  determined  by  means 
of  the  analysis  associated  with  linearized  compressible— flow  theory. 

If,  moreover,  a  Cartesian  coordinate  system  is  used  such  that  the  wing 
is  situated  on  or  in  the  immediate  vicinity  of  the  xy  plane  and  the 
stream  flows  parallel  to  and  in  the  direction  of  the  positive  x  axis, 
it  follows  that  the  basic  equation  for  the  perturbation  potential 
cp(x,y,z)  can  be  written  in  the  form 

P2  ^xx  "<Pyy  -'Pzz  =  0  ^ 

where  02  =  Mo2— 1,  Mo  being  the  free— stream  Mach  number. 

The  application  of  equation  (l)  to  wing  theory  is  essentially  a 
mathematical  problem  Involving  the  solution  of  a  differential  equation 
with  given  boundary  conditions.  Consistent  with  the  assumptions  of 
linearized  compressible  flow  theory,  or  small— pertur bat 1 on  theory,  the 
boundary  conditions  expressing  the  prescribed  physical  conditions  are 
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given  always  at  z  =  0  and,  as  a  consequence,  boundary  conditions  as 
well  as  the  solutions  are  superposable . 

The  techniques  used  in  the  solution  of  wing  problems  are,  for  the 
most  part,  adaptations  of  existing  mathematical  methods  to  the  specific 
type  of  boundary  values  and  their  supporting  surfaces  that  occur  in 
aerodynamics.  In  particular,  it  is  often  possible  in  theory  to  maim  a 
reduction  in  the  number  of  independent  variables  by  virtue  of  known 
geometric  or  physical  conditions.  The  conical-flow-field  analysis  of 
Busemann  (reference  1)  provides  in  this  way  a  means  of  descending  from 
a  three-  to  a  two-dimensional  potential  equation. 

A  conical  flow  field  is  one  in  which  the  perturbation  velocity  com¬ 
ponents  and  the  induced  pressures  are  constant  in. magnitude  along  any 
ray  from  the  apex  of  the  field.  In  this  case,  the  perturbation  poten¬ 
tial  may  be  written  in  the  form 


cp(x,y,z)  =  xfa 


Pyf2(^  — . =  0Zf  (  z  Py  A 
x  x  J  Pjris\  |3y' J3y )  Pzf3  ^  pPp£"  J 


where  cp  is  a  homogeneous  function  of  degree  one  in  the  three  variables 
An  obvious  generalization  of  this  concept  leads  to  the  consideration  of 
homogeneous  potential  fields  of  higher  degree  or,  as  they  are  sometimes 
called,  quasi— conical  fields.  If  cp  is  homogeneous  of  degree  K  +  1 
it  follows  that 

0 

*<W)  -  *(!£)-  <*>■■“  *(,§,£)-  (*>»«  *°(i$)(3) 

Equation  (2)  yields  conical  velocity  fields,  the  degree  of  homogeneity 
being  zero,  while  for  equation  (3)  the  quasi-conical  velocity  fields  are 
homogeneous  and  of  degree  k.  Applications  of  these  quasi-conical 
fields  to  pitching  and  rolling  triangular  wings  have  been  given  by 
Brown  and  Adams  (reference  2),  while  Ribner  (reference  3)  has  used  sim- 
ilar  methods  in  the  consideration  of  cancellation  elements.  Further 
examples  may  be  found  in  the  literature. 

If  new  variables  are  introduced  in  equation  (l)  such  that 


Py  _  _ 
—  - 


=  it  <p(x,y,z)  =  CXK+1  G(q,|;) 


equation  i^8  ^  ar^>^trar^  cons^ant-’  tlle  "transformed  partial  differential 
(b2-i)  %  +  2r)C  %  +  (£2-i)  _  k[2t)  ^  +  &  Og  -  (k+i)  o]»  o  (5) 
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Thus,  for  a  quasi— conical  flow  field  with  apex  at  the  origin  of  the 
coordinate  system,  the  resultant  differential  equation  is  elliptic  for 
all  values  of  satisfying  the  inequality 

q2(;2  -  (1-T)2  )  (l-£2 )  <  0 

that  is,  for  all  points  inside  the  foremost  Mach  cone 

x2  -  P2y2  -  32Z2  =  0 

The  analysis  of  particular  problems  is  therefore  intimately  associated 
with  the  study  of  two-dimensional,  elliptic— type  equations  and  is  espe¬ 
cially  suited  to  the  use  of  complex— variable  theory.  This  is  the 
approach  taken  by  many  investigators.  In  references  4  and  5,  Lagerstrom 
and  Germain  have  developed  these  methods  in  considerable  detail. 

A  different  approach  to  the 
study  of  lifting  surfaces  in  conical 
flow  fields  has  been  given  by  Brown 
(reference  6)  and  in  reference  7. 

In  this  approach  a  basic  lifting  ele¬ 
ment  carrying  a  uniform  load  distri¬ 
bution  and  extending  radially  from 
the  apex  of  the  field  (see  sketch)  is 
considered  first.  The  induced  veloc¬ 
ity  field  is  calculated  for  such  an 
element  lying  in  the  plane  of  the 
wing  and  inclined  to  the  stream  / 

direction  at  an  arbitrary  angle  6. 

The  solution  of  a  particular  problem 
then  proceeds  along  one  of  two  lines. 

If  the  loading  is  given,  the  strength 
of  each  element  is  fixed  and  the  cal¬ 
culation  of  the  lifting— surface  geometry  depends  only  upon  carrying  out 
the  integration.  This  is  referred  to  as  a  direct  problem.  On  the  other 
hand,  if  the  geometry  of  the  wing  is  given,  the  loading  is  unknown  and 
the  strength  of  the  elements  must  be  adjusted  so  that  the  resultant  ver¬ 
tical  induced  velocities  are  consistent  with  the  given  wing  slope  every¬ 
where  on  the  plan  form.  The  solution  of  such  a  problem  depends  upon  the 
inversion  of  a  relatively  simple  singular  integral  equation  and  is 
referred  to  as  an  inverse  problem. 

As  was  pointed  out  in  reference  8,  similar  methods  apply  to  non- 
lifting  problems  in  conical  flow  fields  and  pressure  distributions  cor¬ 
responding  to  conical  elements  of  thickness  can  be  calculated.  In  such 
cases,  however,  the  direct  problem,  that  is,  the  one  involving  the  eval¬ 
uation  of  an  integral,  is  the  one  in  which  the  slope  of  the  wing  surface 
is  given;  and  the  Inverse  problem,  that  is,  the  one  involving  the  solu¬ 
tion  of  an  integral  equation  is  the  one  in  which  the  shape  of  the  pres¬ 
sure  distribution  is  prescribed. 
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Tile  present  paper  is  concerned  with  the  generalization  of  the 
basic  elements  of  references  7  and  8  and  their  application  to  lifting 
and  nonlifting  problems  in  qua si-con ical  fields  of  flow  for  cases 
involving  subsonic— type  leading  edges.  Only  solutions  to  the  inverse 
problems  will  be  considered  and  at  all  times  these  will  be  obtained  by 
inverting  the  integral  equation. 

The  orders  of  singularities  that  arise  in  the  analysis  are  such 
that  it  is  convenient  to  use  the  concepts  of  the  finite  pant  and  gener¬ 
alized  principal  part  of  improper  integrals.  These  generalizations  will 
prove  to  be  of  importance  for  their  notational  efficiency  and  permit  a 
simplified  treatment  of  the  derivatives  of  singular  integrals.  Hadamard 
(reference  9)  has  pointed  out  clearly  the  necessary  steps  in  the  treat¬ 
ment  of  such  improper  integrals,  but  did  not  stress  the  role  of  the  dif¬ 
ferential  operation  in  obtaining  his . integrals.  Since  some  differences 
exist  between  Hadamard* s  definition  of  the  finite  part  and  the  one  used 
here,  when  extensions  to  multiple  integrals  appear,  a  different  notation, 
consistent  with  reference  10  has  been  adopted.  The  generalization  of 
the  principal  part  has  also  been  discussed  in  reference  3. 


LIST  OF  IMPORTANT  SYMBOLS 


a0 

0(d) 

m 


mo 


Mo 

P 

Po 

p 

q. 

Q 


speed  of  sound  in  free  stream 

load  distribution  on  lifting  surface  as  a  function  of  Q 
slope  of  radial  element  relative  to  free— stream  direction 
slope  of  right  leading  edge  relative  to  free— stream  direction 
slope  of  left  leading  edge  relative  to  free— stream  direction 

free— stream  Mach  number 


local  static  pressure 

free— stream  static  pressure 

angular  rate  of  roll  in  radians  per  second 


free— stream  dynamic  pressure 


angular  rate  of  pitch  about  wing  vertex  in  radians  per 
second 
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Vo 

u 

w 

*>7,  z 
a 

0 

& 

So 

&i 

Au 

Aw 


Ap/q 


free— stream  velocity 
streamwise  perturbation  velocity 
perturbation  velocity  normal  to  plane  of  wing 
Cartesian  coordinates  introduced  in  equation  (l) 
wing  angle  of  attack 


angle  between  free— stream  direction  and  line  through  wing 
vertex 

angle  between  right  leading  edge  and  stream  direction 
angle  between  left  leading  edge  and  stream  direction 
discontinuity  in  u  in  plane  of  wing  (uu  —  u^) 
discontinuity  in  w  in  plane  of  wing  (wu  —  w^) 

(p i  -  pu) 

load  coefficient  — — - — 


K 

^  1 
X 

M- 


<P(x,y,z) 


constant  determining  degree  of  homogeneity  of  quasi-conical 
velocity  field 
(See  equation  (3).) 

conical  variables  introduced  in  equation  (^) 
m]3,  hIqP,  miP 

slope  of  wing  surface  relative  to  free— stream  direction 
Mach  angle  (arc  cot  p) 
free— stream  density 

region  of  integration  in  equations  (10)  and  (27) 

perturbation  Telocity  potential  introduced  in  equation  (l) 

function  related  in  equation  (^)  to  perturbation  Telocity 
potential  of  a  quasi-conical  flow  field 
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Subscripts 


u  denotes  conditions  on  upper  surface  of  wing 

f 

l  denotes  conditions  on  lower  surface  of  wing 


ANALYSIS 


This  investigation  is  confined  to  a  consideration  of  inverse  prob¬ 
lems;  that  is,  problems  that  require  the  inversion  of  an  integral  equa¬ 
tion*  As  has  been  pointed  out,  these  problems  correspond  to  the  two 

following  cases:  either  the  load 
distribution  over  a  given  lifting 
surface  is  to  be  determined  or  the 
thickness  distribution  corresponding 
to  a  prescribed  pressure  distribu¬ 
tion  is  to  be  calculated.  The  given 
conditions  must,  of  course,  be  such 
that  a  quasi— conical  flow  results. 
First,  therefore,  the  plan  form  will 
be  chosen,  as  shown  in  the  accom¬ 
panying  sketch,  so  as  to  have  an 
apex  at  the  origin  of  coordinates 
and  to  be  of  semi-inf  ini  t©  extent. 

The  traces  of  the  foremost  Mach  cone 
are  inclined  to  the  positive  x  axis 
at  the  Mach  angle  =  ± arc  cot  0 
and,  since  only  subsonic  leading  edges  are  being  considered,  the  leading 
edges  of  the  plan  form  are  inclined  at  angles  smaller  in  magnitude 
than  }i .  Denoting  these  angles  by  50  and  and  neasuring  them  from 

the  x  axis  positively  in  the  conventional  counterclockwise  direction, 
it  follows  that  the  equations  of  the  leading  edges  are 

y  =  x  tan  =  hiqX  and  y  =  x  tan  Si  =  mix  (6) 

In  the  sketch  SQ  is  positive  while  Si  is  negative. 

The  boundary  conditions  for  the  two  types  of  problems  may  be  stated 
as  follows,  where  subscripts  u  and  l  are  used  to  denote  conditions 
at  z  =  0+  and  z  =  0-,  respectively: 

Lifting  case;  Over  all  the  xy  plane  Aw  =  wu  —  =0  and, 

except  for  the  region  occupied  hy  the  plan  farm,  An  =  Uy  —  u^=0. 

On  the  plan  form,  vertical  induced  velocity  is  specified  in  either 
of  the  forma 
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% 


* 


vu  "  WZ  "  Vo  **  7i  (j)  =  VD  ^  72 

where  7i  or  72  are  known  polynomials  in  x/y* 

Nonlifting  case :  Over  all  the  xy  plane  Au  =  u^  —  u^  =  0 
and,  except  for  the  region  occupied  "by  the  plan  form. 

Aw  =  wu  —  W2  =  0.  On  the  plan  form,  the  streamwise  induced  veloc¬ 
ity  is  specified  in  either  of  the  forms 


«u  -  U1  -  vo  y*  vi  ( |)  =  To  **  v2  (f  ) 

where  v  or  v2  are  known  polynomials  in  x/y. 

The  solution  to  the  two  problems  will  he  attained  after  con¬ 
sidering  first  a  lifting  element  and  a  thickness  element  and  then, 
for  each  of  these  problems,  the  basic  integral  equation  is  deter¬ 
mined  by  summing  the  appropriate  elements.  These  derivations  are 
given  in  the  following  sections.  In  small— perturbation  theory  the 
local  load  in  coefficient  form  is  related  to  u  by  the  expres¬ 
sions 


£p 

1 


Pl~Pu  _  2Au  4uu 


5P0V 


(7) 


and  the  slope  in  the  streamwise  direction  of  an  arbitrary  surface 
z  =  z(x,y)  is  related  to  vertical  induced  velocity  hy  the  expres¬ 
sions 


X 


u 


dzu  =  Wu 

dx  VQ 


dz  i  w. 


> 


J 


(8) 


Lifting  Case 


Upwash  field  of  lifting  element,—  Consider  a  radial  element  emanat¬ 
ing  from  the  origin  and  assume  that  the  load  carried  by  the  element  is 
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(9) 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


/ 


where  C  is  a  constant  for  a  fixed 
position  of  the  element.  If  the 
element  is  inclined  to  the  x  axis 
at  an  angle  6,  its  upwash  field  can 
he  calculated  by  subtracting  the 
induced  fields  of  two  triangular 
plan  farms  with  vertex  angles  equal 
to  6  +  AS  and  5,  each  triangle 
having  one  side  fixed  for  conven¬ 
ience  along  the  x  axis  (see  sketch). 
Assume ,  first,  that  5  is  positive. 

As  shown  in  reference  10,  the  upwash 
field  of  the  triangle  can  be  calcu¬ 
lated  from  the  known  load  distribu¬ 
tion  by  means  of  the  fundamental 
formula 


(x-xi)^  (xi,yi)  dxi 


(y-yi f 


(10) 


where  the  region  t  is  the  area  on  the  plan  form  that  lies  ahead  of 
the  traces  of  the  Mach  forecone  from  the  point  (x,y,0).  The  bars  on  the 
integral  sign  indicate  that  the  generalized  principal  part  of  the  inte¬ 
gral  is  to  be  evaluated  (see  appendix).  By  definition,  if 

=  G(y,yi)  +  constant  (lla) 

(yi-y) 

is  a  known  indefinite  integral,  the  definite  integral  is  evaluated  as 
follows 

f  =  <>(y,b)  -  <*(y,a),  a,b  i  y  (lib) 

In  case  the  singularity  in  the  integrand  lies  outside  the  region  of 
integration,  the  definition  yields,  of  course,  the  conventional  definite 
integral.  In  some  cases  in  the  following  analysis  the  principal-part 
sign  will  be  used  to  express  relations  valid  for  singularities  both 
inside  and  outside  the  range  of  integration. 
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The  integration  with  respect  to  Xi  in  equation  (10)  is  carried 
i  out  over  the  area  bounded  by  the  lines 

7i  =  0,  xi  =  yi  cot  6  =  yi/m,  (i-xi)  =  ±0(y-yi) 

and,  with  the  substitutions, 

0  =  pm,  t|  =  0y/x,  t)i  =  0yi/x  (12) 


leads  to  the  results 


0C 


0(l+T}) 

1+0 


— ^|1 . -gV  (1  -  f  -  (r\-r\if  d.T|ij 

(n^u )  0 


O<0,  — 1  <  TJ  <  s 

(13a) 


0c  r 
knJQ 


Q(l-*l) 

i-0 


ill 


(n-n)‘ 


Ji i- 


ni.s 

e  ’ 


-  (l-1i)  drji;  O<0,  0  < q  <1 

(13b) 


If  0  and  m  are  negative,  the  limits  in  equations  (13a)  and  (13b)  are 
reversed  so  that  as  the  limits  are  now  written  changes  in  sign  are 
»  required  in  the  equations. 

In  order  to  obtain  the  upwash  field  of  the  required  lifting  element, 
it  is  sufficient  to  perform  a  direct  differentiation  for,  if  wu  of  a 
plan  form  with  vertex  angle  &  is  of  the  form 

*u  =  (s.j'n) 

it  follows  that  when  O<0,  the  value  of  upwash  induced  by  the  element 
may  be  denoted  dwu  and  is 

dwu  =  f(0+A9,Tj)  -  f(0,T|)  =  d0 

and  when  0  <  0 


dwu  =  f  (0 ,T)  )  -  f  (0+£0, tj)  =  -  d0 


If  this  process  is  carried  out  in  equations  (13a)  and  (13b)  and  if  the 
trans  format i on 
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is  introduced,  the  expressions  for  dwu  he come 


/pt  dwu  _  pc  de  K  rn 


ft  t  d9  eKf" 

W  v0  iwt  jx 


- - -t-  ^ — .  -  i<  t]<  6 <  i 

(e-ri)(0-t)K+1  t^i-t2 

_ .  -1<9<T|< 


(14a) 


(11*) 


Equations  (l*+a)  axui  (l4b)  provide  the  up-wash  fields  for  any  radial  ele¬ 
ment^  regardless  of  the  sign  of  0,  Integration  by  parts  leads  to  the 
alternative  forms 


&wh 


pc  d0(K+l)  eK  r ^  (r|-t)K  Vl-t2  dt 
*»  J-i  <«-t)K+s  t  ' 


1<  T)<  &<  1 


(15a) 


K 

P  \  dw, 

iy  ~v 


pc  d0 (k+1 )  0K 

tat 


l 


(n-t)K  Vh 


-t2  dt 


/  xK+2 

(0— t)  t 


,  -  i<e<n<i 


(15b) 


0 


* 


where  in  one  of  the  integrals  the  singularity  in  t  requires  the  use  of 
a  Cauchy  principal  part. 

Derivation  and  inversion  of  integral  equation,—  If  now  the  lifting 
elements  cover  the  region  between  0 i  and  0O  and  C  is  a  function 
of  0  determining  the  lift  carried  along  the  radial  element  at  that 
pointy  upwash  produced  by  the  resultant  plan  form  is  given  by  the  rela¬ 
tion 


/'m  wu  _  iim  p(k+i)  [  rn  e  nK  r/nx  .0  r 
K*)v0-^°  s  c(e)a7 


11  (n-t)K  Vi-t2  dt  + 

t  *  x \ K+2  , 


(e-t) 


P  °  eK  c(e)  d9  Pn  —  dt  ] 

Jn+e  J— l  (0— t)K+  t  -I 


TpfG  ^—1  (0— t) 

This  result  can  be  written  as 


(l6a ) 


PiiLti)  r°  e k 

^  ^0i 


C(0)  H(0,tj)  d0 


(l6h) 


4 


u 
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where 


(n-t)K  ^/l-t2 
(0-t)  t 


0 1<!  0<  T) 


(r)— t  )K  f/l-t2 
(0-t  f+£  t 


at. 


n<  eo 


The  function  H(0,q)  has  a  simple  pole  at  Q  -  r\,  and  the  integral 
expression  for  wu  in  equation  (l6b)  is  therefore  evaluated  as  a  Cauchy 
principal  part. 

The  boundary  condition  to  he  satisfied  by  equations  (l6a)  and  (l6b) 
is  that  wu/xK  Yo  a  polynomial  of  degree  k  in  the  variable  tj* 
It  follows  that  the  (k+1)  derivative  of  the  right-hand  member  of  equa¬ 
tion  (l6b)  must  vanish* 


Thus 


0  =  f  °  dK  0(0 h(©,ti)  ae  (17) 

vhere  use  is  made  of  the  generalized  principal  part  of  an  integral 
defined  as  (see  appendix) 

f*3  A(xi)  dxi  _  _1_  / _d  \  pb  A(xi)  dxx 
Ja  (x!-x)n+1  n!  V  W  Ja 


/  ^  n*1+1  b 

fa  A(xi)*nl*-xi|dxi 

Here  again  the  definition  applies  regardless  of  the  value  of  x  but  is 
of  particular  significance  when  x  lies  within  the  region  of  integra¬ 
tion. 


Continuing  the  calculation  in  equation  (17)^  one  has,  after  taking 
the  derivatives  with  respect  to  t], 

0  =  r9°  e*  cle)  aa 

*  *ei  (e-n)K+z 
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which  may  he  written  in  the  form 


v  ^  +  1  9  K 


o  e*  c(e)  cL9 


The  function  C(@)  is  thus  to  he  found  through  the  inversion  of  the 
integral  equation 

K 

r9 o  9K  c(e)  as  r1  , 


The  inversion  of  the  integral  equation 


f(x)  =  r  eLiiii,  a<  x<  b 
Ja  2-1 


is  known  to  be 


g(x)  = 


L=  [a  - 

)  ( T — fl.  )  L  l- 


(h— x)(x-a) 


f(s)y(H)(N) 

i-6 


where  A  is  an  arbitrary  constant  to  he  determined  from  physical  con¬ 
siderations,  Thus,  the  solution  to  equation  (l8)  for  0i<  G<  GQ  is 


9k  C(9)  = 


(0o-0)(0-01) 


n. 


°  i=o 


and  this  leads  to  the  expression 


e  c(e)  = 


v  q—9  )  (d-9 1) 


where  the  coefficients  hi  are  functions  of  9Q  and  9 1  hut  not  func¬ 
tions  of  9 • 


Relation  of  general  solution  to  wing  geometry*—  From  equation  (9) 


the  loading  on  the  plan  form  is,  since  9  =  f3y/x. 
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££  _  (  x 

a  V  0 


K  K+l 


V  b^1 

i=o  ^(^q— 0)(9— ^i) 


where  the  coefficients  bi  must  "be  determined  from  known  information 
at out  the  surface  geometry. 


Consider  next  the  identity 

_ dt _ _  (-1)1 

a  (0— t)i+1*/  (b— t)(t-a)  i! 


L 

,  \  bd  ) 

•  \  / 


a<0<b 


a  (0— t)  v(b-t)(t-a) 


where  i  is  zero  or  a  positive  integer.  This  expression  implies  the 
equality 


f1  ( ^  ^/l— 

f  K+e 

■1  (0-t)  t 


dt  =  0 


since  the  latter  form  can  he  broken,  by  expansion  into  rational  frac¬ 
tions,  into  integrals  like  the  left  member  of  equation  (22).  The  equal¬ 
ity 


^  (t]— t )K */ 1— t2 


-1  (0-t)K+2  t 


(Tj-t)  V  1— t: 


(0-t)K+2  t 


follows  where  the  principal— part  sign  is  needed  on  hut  one  side  of  the 
equation,  depending  on  the  value  of  G  relative  to  r\.  From  this 
result  and  equations  (l6)  the  following  relations  are  obtained: 


0 ) 

wu 

P(*+l) 

x) 

Vo 

4it 

K 

~) 

wu 

3(k+1) 

XJ 

Vo 

k -It 

*  om  *  f  n,  ,<0  (23a, 


K  n 

9k  0(9)  d0  f 
1  ^1 


a  t  (0— t ) 

,T1  (q-t  )K  Vl-t2 
t(0_t)K+2 


dt,  tj>0 


(23b) 


Ik 
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Expressing  the  principal  parts  in 


The  range  of  tj  in  these  two  equa¬ 
tions  has  been  restricted,  respec¬ 
tively,  to  negative  and  positive 
values  in  order  to  avoid  mathemat¬ 
ical  difficulties  arising  when  sin¬ 
gularities  occur  simultaneously  at 
t  =  0  and  t  =  0*  The  shaded  por¬ 
tion  of  the  accompanying  sketch  is 
the  region  of  integration  in  equa¬ 
tion  (23a).  In  the  sketch,  the 
inequalities  -l<0i<O,  O<0Q<^1 
have  been  assumed.  This  implies 
that  the  plan  form  has  two  subsonic- 
type  leading  edges  and  such  a  con¬ 
dition  will  be  assumed  to  apply 
henceforth. 

It  is  particularly  convenient 
to  invert  the  order  of  integration 
in  equations  (23a)  and  (23b). 
the  forms 


71  (ri-t)KvT^  ^  =  lim  1  /  _5^+1  n*  (rj-tf  ^/I~^ 

J-i  t(0-t)K+2  cl->o  (B+i)i  ^  da/  J_±  t(e-t-a) 


£ 


t(0-t)K+2 


dt  =  1Jf 
a— >0 


(r  +1 ) !  \  5a 


K+l 


t(0— t-a) 


substituting  from  equation  (20),  and  inverting  the  order  of  integration, 
leads  one  to  the  forms 


lim  — ft  ( _d_  \  (n-t  f  <J_  1— t2 

a^°  k*Kl  \  5a  )  J_x  t 


K+l 

I 


i=o 


bi  0i  d9 


(0— t-a)  >J(q0-q)(q-Q 


TJ  <  0 


4 


* 
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3  )  =  11m _ 3_  f  jL  1 

x/  V0  a->°  k*K i  \&tj 


K+l  1 


«  /TI5^  nd 


p  l  (n-t)  vi— t2dt 


y  Tij  y  8^  j(t-ct,)j~ide 

°  1=0  M 

*/ (do~& )  (0— 0i ) 


'®°~a  (tj— t)K  </l— t2dt  f^o 


]T  bi  (t+a) 


0i  (0-t-a) J (0Q — 0 ) (0 — 0 1 ) 


^  bj  (t+a)1  d0 


13  (r]— t)  Vl— t2dt  l=o  _ 

So-00  t  (0 — t — a)  J(e0-e)(e-ei) 


The  double  Integrals  occurring  first  in  the  rigiht-hand  members  of  these 
equations  are  of  degree  K  in  a  and  their  derivatives  consequently 
vanish*  Moreoever.,  from  the  identities 


a  (0-x)  */ (b— 0 )  (0— a ) 


(x-a. )  (x— b ) 


(b— x)(a— x) 


j  b<  x 


, a<x<b 


,  x<a 


it  follows  that  the  last  terms  in  the  right  members  vanish  and  that  the 


expressions  for  wu  become 
/  \K  f 

(  3  )  25i  =  lim  _±(  ^ 

\  VQ  a->o  4k  I  \  aa> 


l9l-“  („-t 

t 


bj(t+a) 


i=o  v  (0q— t-a)(01-t-a) 


dt.,  tj  <  0 
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x  nk  ,  x  U+i 

/P]  =  11m  _J_ 

\xj  yo  a~> o  4k J  \  da ) 

9c-  (T)_t)Kv/H? 

t 

x 

If,  as  in  the  appendix,  the  finite  part  of  an  integral  is 

Tb  A(xi)dxi  =  21  A(xi)dxj.  = 

(x1-b)i+l/2  1.3...(2i— 1)  '-'a  'db'  (xi-h)1/2 

gl  'j  f°  ^2^1 

1.3... (21-1)  \db'  ^a  (xi-b)1/2 
the  expressions  for  wu  may  he  written  in  the  form 


K+l 


I 


b^(t+a)^ 


l^Q  ^(t+a^0)(t+a-91) 


dt. 


rj  >0 


K+l 


I V1 

1=0 

-T1—  It,  T)  <0 

v  (eo-t)(0i-t) 


(25a) 


(  P  ^  =  _P_  /0°  (Ti-t)Vl-t2  / _d_\" 

\ x/  VQ  4k  J  t  \at/ 


K+l 


I 

i  =n 


hjt1 


*/  (t— 0O )  (t— 01 ) 


dt,  tj>0  (25b) 


Equations  (25a)  and  (25b)  are  the  fundamental  equations  for  a 
lll'H^S  surface  and,  since  wu/^o  4ias  been  assumed  known  as  a  polyno¬ 
mial  in  q  =  Py/x,  it  remains  merely  to  determine  the  unknown  coeffi¬ 
cients  bj  by  equating  coefficients  of  rj  on  both  sides  of  the 
expressions.  In  the  form  given,  it  appears  that  for  K>0  the  number 
of  equations  obtainable  exceeds  the  number  of  undetermined  coefficients. 
No  general  theorems  of  determinancy  have  as  yet  been  established  as  to 
the  uniqueness  of  the  solutions,  but  applications  to  be  made  later  will 
indicate  the  techniques  involved  in  calculating  specific  examples.  In 
the  laterally  symmetrical  case,  where  — 0i=0o,  solutions  are  easier  to 
determine  and  the  fundamental  equations  are 
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_±  r ~^°  (n-t)K 
itK  !  J_x  t 


K+l 

V 


1=0 


dt,  t]  <0 


p  r9°  (ti-t  f  Vi-t2 

^KiJi  t 


i=o 


at,  t)  >0 


(26a) 


(261) 


Nonlifting  Case 


A  radial  element  emanating  from  the  origin  is  to  he  constructed 
such  that  it  has  a  quasi-conical  thickness  distribution 

V  =  c/ 

where  C  is  a  constant  and  is  the  streamwise  slope  of  the  element 

as  defined  in  equations  (8).  The  derivation  of  the  induced  pressure 
field  associated  with  the  element  follows  closely  the  analysis  in  the 
lifting  case.  Thus,,  a  triangular  plan  form  is  first  considered  where 
one  side  is  parallel  to  the  stream  direction  and  with  a  vertex  angle  6. 
From  reference  10.,  pressure  coefficient  can  be  written  in  the  form 


C 


P 


_ (x-xi)  dxi _ 

[(x-Xx)2H32(y-yi)2]3/2 


(27) 


where  the  region  t  is  the  area  on  the  plan  form  that  lies  ahead  of  ‘the 
traces  of  the  forecone  from  the  point  (x,  y,  0)  and  the  integration 
with  respect  to  xi  involves  the  finite  part. 

In  the  notation  of  equations  (13)*  the  analogues  to  equations 
(15a)  and  (15b)  are 


2  c  ae  9k  (k+i  )  t(ii-t)K  at 

J-i  (e-t)K+zJi~^ 


ac  =  -  2  G  d9  0k(k+1) 
P  itp 


t(n-t)  at 
(e-t  )K+2  Vi-t2 


— 1<  T)<  e<  1 


(28a) 


X 


3 


-1<0<TJ<1 


(28b) 


NACA  TN  2^97 


19 


If  C  is  a  function  of  0  and  the  thickness  elements  cover  the  region 
"between  9 i  and  0O,  pressure  coefficient  on  the  plan  form  is  given  by 
the  expression 


0)‘ 


_  lim  ~£(K+l) 


P  €->o 


eKc(e)do  ^  -  + 

0i  J\  (6— t)K+2'/l-t2 


0KC(0)d0 


t (tj — t)K  dt 
^-i  ( 0-t  )K  +2  Vl— t2 


The  "boundary  conditions  require  pKC_/xK  to  be  a  polynomial  of 
degree  k  in  tj.  If  the  (n+l)st  derivative  of  equation  (29)  is  se" 
equal  to  zero,  the  relation 


0  =  0KC(0)d0  ( 1 

Jdi  \  <*1 


t(r]-t)K  dt 
!  (0-t)K+2^t2 


'°  a«c«,)da  (>T  f*  W** 

\^T)/  J_1  (9— t  )K+2  «/ 1— t2 


holds,  and  after  further  differentiation  reduces  to 


qkc(o) da  /  a 


0i  (0— TJ  ) 


l  r  o  e*C(e)d6> 
(K+l  )]  J Q±  (0—^  ) 


The  function  C(0)  satisfies  the  same  integral  equation  that  arose  in 
the  lifting  case.  (See  equation  (l8).)  The  solution  can  therefore  "be 
written,  as  in  equation  (21),  in  the  form 

.  vK  k+i 


V  -  o(e)  7"  -  (|)  X 


i=o  V(9q-^)(9-9  1) 


The  equivalence 


f11  ~b(n— t ) K  dt  =  t  (q— t  )K  dt 

-1  (a_t)K+2-/i-t2  Ji  (e-t)K+2  Vi-t2 


(30) 
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permits  the  rewriting  of  equation  (29)  in  the  forms 


(31a) 

(31b) 


Substitution  for  0KC(e)  from  equation  (30),  inversion  of  the  order  of 
integration,  and  use  of  equations  (24)  leads  to  the  fundamental  rela¬ 
tions 


K+l 


K+l 


(t— 0O) (t— 9i) 


dt,  t]  >0 


When  =  -B0,  these  equations  become 


(32a) 


(32b) 


The  determination  of  the  thickness  distribution  corresponding  to  a 
given  pressure  distribution  can  thus  be  obtained  from  the  above  equa¬ 
tions  by  equating  coefficients  of  t\  and  solving  for  the  unknown 
coefficients  a*j_.  Specific  examples  will  serve  to  make  the  steps 
clearer;  such  problems  will  be  considered  in  the  following  section. 
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APPLICATIONS 


Since  homogeneous  fields  of  low  degree  have  already  received  con¬ 
siderable  attention,  several  results  have  been  published  previously.  In 
the  case  of  thickness  problems  corresponding  to  specified  pressure  dis¬ 
tributions,  however,  solutions  have  never,  so  far  as  is  known,  been 
sought  in  terms  of  the  given  pressures.  Rather,  the  thickness  has  been 
assumed  known  and  the  resulting  pressure  distribution  calculated.  This 
latter  attack  involves  no  question  as  to  uniqueness  and  a  one-to-one 
correspondence  certainly  exists.  When  pressure  is  prescribed  first, 
however,  it  becomes  necessary  to  consider  the  possibility  of  nonunique— 
ness.  In  two-dimensional,  low-speed  flow  a  freedom  of  choice  is  known 
to  exist  and  leads  to  the  introduction  of  purely  circulatory  flow  which, 
in  turn,  provides  the  mechanism  of  lift.  No  analogue  to  this  occurs  in 
the  low— speed,  two-dimensional,  nonlifting  case  when  the  body  is  smooth 
and  is  assumed  to  close.  In  the  following  developments  a  multiplicity 
of  solutions  will,  however,  occur  in  the  nonlifting  case  and  bodies  with 
given  pressure  distributions  retain  a  degree  of  freedom. 


Equations  (25)  and  (26),  together  with  equation  (21),  suffice  for 
the  solution  of  qua si-coni cal  lifting  problems  while  equations  (32) 
and  (33 )*  together  with  equation  (30),  apply  to  symmetrical  wings.  In 
the  applications  to  follow,  the  division  into  lifting  and  nonlifting 
cases  has  been  maintained.  The  detailed  treatment  of  equations  (26)  and 
(33)  can  be  further  simplified  if  the  problems  are  separated  into  cases 
involving  symmetry  and  antisymmetry  about  the  x  axis  of  the  imposed 
boundary  conditions.  Suppose,  first,  that  the  given  values  of  upwash 
and  pressure  coefficient  in  these  equations  are  odd  functions  of  q. 

It  follows  from  physical  considerations  that  the  loading  or  surface 
slope,  respectively,  will  be  an  odd  function  of  q  and  that  conse¬ 
quently  the  unknown  coefficients  ai  or  must  vanish  for  even 

values  of  the  subscript  i.  If  the  transformation  t  =  — r  is  made  in 
either  equation  (26b)  or  (33b),  the  pairs  of  equations  (26a)  and  (26b) 
or  (33a)  and  (331)  yield  consistent  sets  of  simultaneous  linear  equations 
that  can  be  obtained  from  single  equalities.  Hence,  for  wu  or  Cp 
expressed  in  odd  powers  of  q. 


+0  £dO  (t|— t)K*/  1— tg 

4kJ  J[  t 


[K/2] 

y  k  ,  t2 

Z_j  2J+1 
J=Q _ 

'/t2-0o2 


(3*0 


or 
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where  the  notation  [*/2]  in  the  summation  denotes  the  largest  integer 
c  ont a  ine  d  in  k/2  . 


In  the  same  fashion,  a  simplification  can  be  achieved  for  symmetri¬ 
cal  boundary  conditions.  Again  using  the  bracket  notation  to  indicate 
the  largest  integer,  the  resulting  equations  become,  when  wu  or  Cp 
are  expressed  in  even  powers  of  tj. 


+3  P9o  (ij-t  )K  v4-t2 
4k  I  t 


[  (k+i  )  /2  ] 

i  ^ 

J=o  ■ 


t2<3 


dt 


or 


g  r9o  t  (ti~ t  )K 


[ (k+i )/s] 


J=o 


t£J 


dt 


(36) 


(37) 


Wings  With  Load  Distributions 


The  vawed  triangular  wing.—  Equations  (21)  and  (25)  lead  directly 
to  the  determination  of  angle-of-attack  loading  on  a  yawed  triangular 
wing.  This  solution  is  well  known  and  was  calculated  in  reference  7  by 
a  method  which  was  a  particular  case  of  the  present  theory. 

The  boundary  conditions  are  that  wu  =  —  VQa  on  the  plan  form, 
hence  k  =  0  and,  from  equation  (21), 

£p  =  bp+bxS 

^  J (9  0-Q  )  (OS i ) 
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Equations  (25)  lead  to  the  equalities 


p  Pdo  Vl-t2  S  b0+bjt 
b-Ji  t  St  </(t-e0)(t-e1) 


_  a  =  z£  £dX  - 1-1:2  A  ^0+^  dt 
4  t  St  ^(0o_t)(0i_t) 

After  integration  by  parts,  these  relations  become 


0  pdo  _ (bp+bjt )  dt 

^  J  i  t2  J  ( 1— t2  )  ( t  0  o )  ( t— 0  !  ) 


/  1  (bo+bjt )  dt 

— 1  t2v'(l-t2)(0o-t)(0a-t) 


Hence 9 


where 


a  -  7"  lo(0oj  ^i)  +  ^1(^0*  ®i 

4  _ 


a  =  4  “~bi  ki(— 9i*  -*o) 


LO(0O,  ^l)  ~ 


t2  v  (l— t2  )  ( t— 90 )  (t— Oi ) 


^i)  — 


4  t  a/ ( l—t2  )  ( t— 0  ©  )  ( t— 0 1 ) 


It  is  apparent  that  equations  (39)  can  he  solved  for  hQ  and  hi 
in  terms  of  the  functions  Lo(0o*  9i)  and  Li(0o,  Si)*  Substitution  of 
the  values  into  equation  (38)  leads  to  the  expression  for  load  coeffi¬ 
cient 


=  2a  /  2G  r  (9q+9i)  9-29q0i 
q  0E  J  Q0-Q\  L  ^(0  -0)(0-01) 
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vhere  E  is  the  complete  elliptic  integral  of  the  second  kind  with 
modulus  Vl-G2  and 

q.  _  I—QqQ i  —  V(l-9og )  (l— 9 12 )  _  C^/(l~ Qq ) (l+8i )  —  */(l— 9i ) (l+9p )  J  ( ) 
9o~ei  2(eo-0!) 

In  the  particular  case  when  =  — 9OJ  the  value  of  G  heconBS  90 
and  the  resultant  loading  on  the  unyawed  wing  be  cobbs 

£p  ^  ka  9qz 

4  PE  V902-92 

where  the  modulus  of  E  is  «/l-eo2.  This  latter  problem  could,  of 
course,  have  been  solved  directly  from  equation  (36)  in  a  much  simpler 
fashion. 

The  rolling  triangular  wing.-  Consider  next  the  case  of  an  unyawed 
triangular  wing  rolling  about  its  axis  of  symmetry.  If  the  angular  rate 
of  roll  is  P  radians  per  second,  the  boundary  conditions  on  the  wing 
are  that  wu  =  —  Py.  In  this  case  k  =  1  and,  from  equation  (21),  the 
loading  on  the  wing  is  given  by  the  expression 


^2.  _  x  bi9 


(42) 


where  the  coefficients  bQ  and  b2  can  be  deleted  since  the  loading 
must  obviously  be  antisymmetrical.  Since  [k/s]  is  equal  to  zero  in 
equation  (34),  bi  is  found  to  satisfy  the  relations 


_L  =  H3bi  fo  f  _d  ^  t 

v0  4  Jx  t  U  /  yt2-o02 


(43a) 


0  =  ba 


(43b) 


Performing  one  differentiation  with  respect  to  t,  in  each  of  the 
two  integrals,  and  then  integrating  by  parts,  leads  one  to  the  expres¬ 
sions 


Vo 


Pti9o2  f9o 
4  Jx 


_ _ _  ffiiQo  d  r  _ dt _ 

t2Vi-t2  (t2-9Q2)3/2  4  t2v/(i-t2)(t2-0o2) 
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r9o 


0  -  1D10O2  T 


t  dt  v  „  d 

A  (t2-S  2)3/2 


f 


t  dt 


/(l-ts)(t2-S  2) 


O  t  — u  J.  "v  \-*-  «  /  \  w  -^o 

The  last  equation  reduces  to  an  identity  while  the  former  one  becomes 


p 

-Pbi9o 

3  ( 

6  E  > 

,  _  PDi 

'  2E-0O2  (E+E)  ■ 

II 

u° 

4 

*o' 

il 

S' 

°  IV) 

L  d^o2)  -1 

(44) 


where  the  modulus  of  the  complete  elliptic  integrals  E  and  E  is 
k  =  VI-Qq^.  From  equations  (42)  and  (44)  the  loading  is 


<1 


4p  0o 2  X  9 


p2Vn  (  E - ^2- 


!^o2 


l^o2 


K 


)  A 


(45) 


2-^2 


and  is  in  agreement  with  the  results  of  reference  2. 


The  pitching  triangular  wing.—  If  an  unyawed  triangular  wing  is 
pitching  about  its  vertex,  the  boundary  conditions  become  vu  =  — Qx 
where  Q  is  the  rate  of  pitch  in  radians  per  second.  From  equa¬ 
tions  (21)  and  (36),  the  loading  and  the  relation  involving  the  undeter¬ 
mined  constants  are 

Ap  _  x  'b0+'bs92 
q  P  v'Qq2^ 

and 

— PQ  _  P  r9°  (t»-t)  ^l-t2  /  a  N*  b0+b2t2  ^ 

Vc  4  ^  t  \dt  /  yts-0o2 

since  K  =  1  and  the  loading  must  be  symmetrical  about  the  x 
Equating  coefficients  of  r\  in  the  latter  equation,  one  has 

axis. 

a  _  r9°  ^i-t2  f  a  f  bo+bst2  dt 

^1  t  \  dt/  <s/t2-0o2 

(46a) 

(  h  1  4t 
v0  b J x  V at/  vt2-0o5 

(46b) 
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or,  after  carrying  out  one  differentiation  and  integrating  by  partB, 


0  =  b2 

4q 


a  r9° 
&o  J 1 


t  dt 


yri-t2)(ts-0O2) 


-  (bo+2b2e02) 


=  — b2 


—  A 


t4  dt 


e.  JU^W^T) 


(bo+2b20o2 )  fX 


t2  dt 


V(i-t2)(t2-a02) 


a  p9° _ dt _ 

^O  Jl  t  </  (1-t2 )  (t2-0o2  ) 


The  first  term  in  the  first  equation  is  zero  and,  after  integration  of 
the  terms  in  the  second  equation,  the  following  relations  are  obtained 


t 


0  =  V  20o2  b2 

=  b0(E^)  +  (E-£02K)b2 


where  the  modulus  of  the  complete  elliptic  integrals  is  k  =  *f\. -0O2 . 
The  solutions  of  these  simultaneous  equations  are 

+8  B£\? Q  ^  -4  k2^ 


*o  = 


V0  [  0o2K+(l-20o2)E  ] 


,  t>2  = 


[  0  02K+  ( 1-20  Q2  )E  ]  V, 


and  the  resultant  loading  on  the  wing  is 


1 


4x  Q 


2  0O2  -  e2 


2-02 


( 


-£oL  K  +  iz££oiE 

l-^o2  l-^o2 


) 


(47) 


(48) 


Differentially  deflected  triangular  wing.-  If  the  two  sides  of  an 
unyawed  triangular  wing  are  deflected  differentially,  vertical  induced 
velocity  on  the  wing  is 


NACA  TN  2497 


27 


and  loading  must  therefore  he  asymmetrical.  It  follows  that  K  =  0 
and,  from  equations  (21)  and  (34), 


Ap  _  hi  0 

<1  Je02-e2 


and 


—a 


=  p  re° 
-4  j; 


Jim 


_a_ 

St 


hi  t 


dt 


An  integration  hy  parts,  in  the  latter  expression,  reduces  the  integral 
to  a  standard  form.  A  further  integration  leads  to  the  equality 


hi 


8  0O a 


and  load  distribution  is  therefore 


Ap  8a  0O  0 

q  0it  Je0£-ez 


(49) 


Triangular  wing  with  parabolic  twist.—  Consider,  finally,  an 
unyawed  triangular  wing  twisted  symmetrically  such  that  its  vertical 
induced  velocity  is  of  the  form 


2u 

Vo 


r  y2 


where  r  is  a  fixed  constant* 


Since  k  =  2*  the  relations 


f X  ^  ho+bg9g 


and 


£  p9o  (n-t  )2  Jl-t2  /  J)f  hp+hgt2 
8v7i  t  '  St/  Vt2-0O2 


apply*  Detailed  analysis  will  he  omitted  in  this  case  since  it  follows 
the  same  pattern  of  development  used  in  the  earlier  cases*  Three  simul¬ 
taneous  equations  involving  the  two  unknowns  hQ  and  h2  are  obtained* 


28 


NACA  TN  2497 


tut  the  equation  relating  coefficients  of  tj  can  he  shown  to  vanish 
identically  and  the  two  remaining  equations  yield  a  unique  answer.  The 
expression  for  load  distribution  is 


AP 

<1 


4  [290gK-(l+90g  )e]+[-(3-902  )K+4(2-0o2  )eJq2 

\_-590*X?+890z(l+eQs  )KE+(49O4-190O2+4)E2]  Jdz-Qz 


(50) 


where  the  modulus  of  the  complete  elliptic  integrals  E  and  E  is 

Vl-90  • 


Wings  With  Thickness  Distributions 


Triangular  wing  with  uniform  pressure.-  In  reference  11,  Squire 
considered  certain  thickness  distributions  for  symmetrical  nonlifting 
wings  in  conical  flow  fields  and  calculated  the  resultant  pressure  dis¬ 
tribution.  The  first  of  Squire's  examples  was  a  triangular  plan  form 
with  a  uniform  pressure  distribution.  It  is  instructive  to  consider  the 
inverse  of  this  problem  and  to  seek  the  wing  assuming  the  pressure 
variation  known.  The  plan  form  is  symmetrically  disposed  with  respect 
to  the  x  axis  or  stream  direction  while  the  boundary  conditions 
require  that  Cp  is  a  constant  over  the  entire  wing.  In  this  case, 

K  *  0  and  from  equations  (30)  and  (37)  the  following  relations  hold 


<*zu  _  ^  _  ap 
41  ^  */0o2— 02 


,  _  2  re°  t  a  aQ 
p  “  p  £ 


(51) 

(52) 


The  value  of  aQ  may  be  found  in  a  manner  quite  similar  to  the  one 
used  in  the  previous  examples.  Thus,  after  the  differentiation  is  per¬ 
formed,  equation  (52)  becomes 


_ t2  dt 

(t2^027^ 


ap  _a_  r9°  t2  dt 

9o  *9oi 


PCp 

2 


^(IIt^)(t2V) 
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After  integration,  this  yields 


Pcp  _  ao  d  E  _  ap(K-E) 
2  9q  beQ  i-eQs 


(53) 


vhere  the  modulus  of  the  elliptic  integrals  is  vl-0o2.  If  aQ  is 
eliminated  from  equations  (51)  and  (53)*  the  slope  of  the  ving  on  the 
upper  surface  is 


azu  =  gcp(i-e02)  i 

2 (K-E )  «/a02-92 


Cp (l—Qp2 )  z 


2(K-S) 

The  ordinate  of  the  upper  surface  results  from  the  integration 


and  is 


dz 


u 


dx 


dx 


Ju 


=  l1!9?2  )c£  =  ■^r9?2^*  Jqq2^2 

2eP(E-£)  20m2  (Z-E) 


(54) 


Triangular  wing  vlth  linear  pressure  gradient*—  It  is  now  proposed 
to  determine  the  thickness  distribution  for  an  unyawed  triangular  wing 
for  which  pressure  varies  linearly  in  the  streamwise  direction.  Set¬ 
ting  Cp  equal  to  bz,  it  follows  from  equation  (30)  that  the  slope  of 
the  upper  surface  is  ezpressihle  as 


z  ao+a202 

a*  P  J @02— s2 


(55) 


since  the  solution  is  obviously  symmetrical  about  the  stream  azis. 
Since  k  =  1,  equation  (37)  becomes 


a0+aat2 


dt 


(56) 


Equation  (56)  is  an  identity  and,  after  the  coefficients  of  t)  are 
equated,  leads  to  the  two  relations 
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0  =  S9°  t  (  3  ^  ao^2,t2  dt 

Ji  /Pt2  \st/  vt2-^02 

(57a) 

02b  _  _r9°  t2  f  S  \  a0+a2t2 

2  X  Vl-t2  Voty  Jt2-4QZ 

(57b) 

It  is  of  interest  to  compare  the  series  of  equations  Just 
developed  with  the  corresponding  equations  in  the  problem  of  the  pitch¬ 
ing  triangular  wing,  Formally,  the  algebraic  steps  are  the  same  and  it 
is  to  be  expected  that.  Just  as  in  the  case  of  equations  (46),  tvo 
simultaneous  equations  will  be  obtained  and  that  their  solutions  will 
provide  the  constants  ap  and  a2*  In  the  present  case,  however,  equa¬ 
tion  (57a)  can  be  shown  to  vanish  identically  and  as  a  result  only  one 
equation  in  tvo  unknowns  remains.  This  means  that  an  infinite  number  of 
possible  solutions  exists.  The  following  calculations  will  supply  the 
necessary  details  to  confirm  these  remarks. 


Consider  equation  (57a)  and  introduce  the  transformation 


t2  =  T,  0O2  =  T  Q 


Then,  by  means  of  the  relations 


St  St  St2  \  Si2  St  / 


the  equation  1)6001168 
°-aoJT 

a2  X 


dT 

f  2t  t  + 

M  1  + 

\  Sr* 

St  J  ,/t-Tq 

dT  ( 

S2 

S  \  T 

VPF  ' 

\  St2 

St  J  */t_t 

=  a, 


r  *2  , 

rT0  2T  dT  s 

pr°  dT  1 

L  std2  J 

!  y(l-T)(r_T0)  Sro  j 

1  J  (1-T  )  (  T-T0  )  _ 

a2 


2t, 


'  X  r°  2T  JHsir-A.  rT°/^4T. 

St0  J1  V  l — t  St  0J1  */  1— t 

d2  pT° _ t  dT _ t  S  rro  dT 

^To  J (1-t)(t-t0)  °  Sto  'J 1  •J  (l— t)(t— 


To)  J 


(58) 

(59) 
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Each  of  the  above  integrals  can  be  evaluated  directly  and  the  bracketed 
terms  are  in  both  cases  zero.  Since  similar  integrals  occur  in  problems 
o£  this  type,  however,  it  is  worthwhile  to  give  the  following  general 
formulas  (n  an  integer) 


-====J===z  =  n  =  O  j  =  -  i(l+To),  n=l; 

l  V(1-t)(t-t0)  2 


=  __rt_  (2n-lH (Tn+l) 

2zn~1  L  (n-l)ini 


.1 — X 

Y  — 2  * —  T0n_J  ,  n  >1  (60a) 

Li  (n— j )  I  j  *  ( J—l ) l (n-J— 1 ) 1  J 


dr  =  —  Jt,  n  =  0  ;  =  |(  T02  +  2T0_3  ),  n=l; 


it(r0— l)  f  (2n-l)V  .  ( 2n+l ) ! 


22n(n+l ) 


(2n-l)Ir 
-  (n-l)inl 


2  n!  n! 


1  _  1 

\  (2 j+1  )!(2n-2 j— 1  )i  T  n— J  J  n>l  (60b) 

L  (n— j)l jljl (n— J— l)i  °  J 


It  remains  to  calculate  the  terms  in  equation  (57b).  If  the  dif¬ 
ferential  relations  in  equation  (59)  are  used,  the  desired  expression  is 


-e2b 


(25^Bl  + aTo®2) 


a2  (  2  i-gDa-  5  JLDi  +  D2 
\  Sto 
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where 


=  f  °— 7  T  dT  =  -  2  I*  dn4  udu  =  -  T  2g  _  1  (l+T  2)E 

J1  a/t(1-t)(t-t0)  Jo  3  °  3  ° 

O  aK 

=  /  —7=  l-dI.  ==  =  -5  /  dn2  udu  =  -2E 

J1  t)(t— tq)  Jo 


d3  =  rT°  t3  ^ 

Ji  ^/t(1-t)(t_to; 


dn2  udu  =  -2B 


dn6  udu  =  2.  r0  ife  +  ^  (l+r0)  Di 


and  the  modulus  of  the  elliptic  integrals  E  and  K  is  */i-Tq* 
Direct  calculation  gives  for  the  coefficient  of  a0  the  expression 

B  (3— To)  -  2K 
(1-r  of 

and  for  the  coefficient  of  a2 

2  t0  E  -  t0  g  -  T02  g 

(1-T0f 

The  value  of  h  from  equation  (6l)  thus  establishes  for  pressure 
coefficient  the  value 


C  =  g:£ 

p  02(i-eo2)2 


ao  2g  —  E(3-0O2 )  + 


ae  [  (0O4+0O2 )  E  -  20 o2  E  J 


where  the  modulus  of  the  elliptic  integrals  is  */l-QQz. 

From,  equation  (55) j  the  ordinate  of  the  upper  surface  of  the  wing 


dzu  1  px  ao^  +a202y2 


202mo2 


/m02r2-y2  +  -  (ao+202mo2a2 )  arc  cosh  22E 

2mo30  7 


(63) 
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In  reference  11,  Squire  considered  the  thickness  distribution  that  is 
obtained  by  neglecting  the  arc  hyperbolic  function  in  equation  (63). 
His  results  correspond  to  the  case  when  a0  is  -2p2mo2a2  and  are 
specifically 


apx 

202mo2 


a/m 


(64a) 


C 


P 


aox 

p2(l-602)2 


[(3 -90*)  K  -  (4— 20o2)  E] 


(64b) 


If  the  wing  is  cut  normal  to  the 
stream  direction  to  form  a  trailing 
edge,  a  triangular  plan  form  and  an 
elliptical  cross-section  result  as 
shown  in  the  accompanying  sketch* 

If  the  root  chord  of  the  wing  is 
cc  and  the  maximum  thickness  at  the 
trailing  edge  is  t,  the  constant 
aQ  in  equation  (64a)  is  equal  to 
02mot/cc^  and  the  analytical 
expressions  for  the  upper  surface  of 
the  wing  and  pressure  coefficient 
are 


\ 


It  is  apparent  that  a  multiplicity  of  thickness  distributions  with 
the  same  pressure  distributions  must  exist.  Consider  first  the  case 
when  aQ  in  equations  (62)  and  (63)  is  zero.  The  surface  ordinate  and 
pressure  coefficient  are,  respectively. 


arc  cosh 


fflpX 

y 


(66a) 
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„  2&2X  e02  r. 


(66b) 


Sr  ?r:;on  of+3w\ving  is  shown  far  * ***■<****  values  <* 

100  ,  f  n  nex^  sketch  (the  curve  denoted  hy  n  =  »)*  Alone  th#=> 

‘te  i.  zero  while  the  »!>.  tblL.™  po£SL 

occurs  at  the  value  of  Dq x/y  satisfying  equation 


mox/y  .  mrc 

— —  1  '  —  =  arc  cosh  ~ 

2  y 


or  at  approximately 


“epe/y  -  1»3128  or  y/mox  =  O.7617 

tloalTr^nL”^48  °f  '4mtIOnS  <65b)  ^  ^  ”  ““ 


as  =  _mot  [(3-9o2)  E  -  (jj-gQn2)El 
P2  20o2co2  [(i+0o2) 


holds.  From  equations  (65a)  and  (66a)  it  follows  that  the  surfaces 


yCTT 

\cn/  V  \nux 


=  /  7  j 

t  ^nioCo/ 


t(3-«02)  E-  (1*-£0o2)e]  nux 

771  p  ,  - : - arc  cosh  -2- 

1(1+0  2 )  E  —  2E] 


(67a) 


(67b) 


Pr°BSUr9  41strl'’-‘tl°“  “  *»■  also  any  ewrf.ee 
2Z  2Zj.  +  2nZg 

T  =  t (l+n )  (67c) 

where  n  is  an  arbitrary  multiplicative  factor. 
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Cross  sections  in  the  x  =  c0  plane  of  surfaces  given  by 
equations  (67a)  and  (67b)  are  shown  in  the  accompanying  sketch  for 


Mq  =i/iT  and  mo  =  l/2.  Also  included  are  sections  calculated  from 
equation  (67c)  for  n  =  1  and  n  =  l/5. 

Since  bodies  with  the  same  pressure  distribution  can  be  found,  the 
possibility  of  combining  results  and  getting  a  body  inducing  no  change 
in  the  free— stream  pressure  should  be  investigated.  From  equation  (62), 
it  follows  that  for  Cp  =  0,  the  arbitrary  constants  aD  and  must 

satisfy  the  relation 

ae  =  _  2K  -  E  (3-0p2 ) 
aD  90z[(l+90s)  K  -  2E] 

and,  from  equation  (63),  the  ordinate  of  the  resulting  surface  is 
expressible  in  the  form 


where 

P3  =  P2/P1 
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and 

Pi  =  (1+9 02 )  K  —  2E 

Ps  =  (3-©02)  K  ~  (^o2)  E 

The  surface  produces  a  real  wing  only  if  2^  remains  positive  and  this 
condition  leads  to  the  inequality 

(m 0x/j)  yCrngx/y)2  -  1  /  arc  cosh  m^/y  ^  P3 


The  range  of  diqZ /y  is  from  1  to  »  and  it  is  easy  to  show  that  the 
left  member  has  the  lower  limit  1  at  niox/y  equal  to  1.  On  the 
other  hand,  P3  can  be  written  in  the  form 

p3 - * - 

P2-2(1-902)(K -E) 

and,  since  K  —  E>0  for  90<1,  it  appears  likely  that  P3  is 
greater  than  one.  A  more  detailed  check  shows,  in  fact,  that  P3  lies 
between  3  and  9  from  which  it  follows  that  the  inequality  can  never 
be  satisfied  in  the  neighborhood  of  the  leading  edge  and  no  real  wing 
with  zero  pressure  coefficient  is  possible.  The  variation  of  Pi,  P2, 
and  P3  is  shown  in  the  accompanying  sketch  in  which  the  variables  are 
plotted  as  functions  of  9C2. 


9  6  1.2 


P  P  P 
3  e  / 

7  4  .6 


5  2  4 


3  0  0 

0  .2  .4  gg  .6  .8  1.0 
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CONCLUDING-  REMARKS 


It  has  been  shown  that  the  assumption  of  quasi— conical  flow  in  a 
supersonic  field  transforms  the  basic  partial  differential  equation  for 
the  perturbation  potential  to  an  elliptic-type  equation  in  two  Inde¬ 
pendent  variables  throughout  the  region  inside  the  Mach  cone.  It  is 
therefore  not  surprising  that  solutions  of  wing  problems,  for  both  the 
lifting  and  the  nonlifting  case,  lead  directly  to  the  consideration  of 
an  integral  equation  (equation  (l8))  of  the  type  encountered  in  two- 
dimensional  subsonic  theory.  In  the  analysis  and  the  applications  of 
this  report,  it  is  shown  that  for  a  large  class  of  specified  conditions 
the  known  inversion  of  the  integral  equation  produces  solutions  that 
require  straightforward  integrations  and  the  solving  of  simultaneous 
linear  equations. 

An  unusual  feature  of  the  resulting  theory  is  the  fact  that  a 
multiplicity  of  solutions  may  appear  in  a  given  problem.  In  retrospect, 
this  degree  of  freedom  is  not  surprising  since  it  is  well  known  that  a 
null  solution  exists  In  two-dimensional  subsonic  theory  and  appears  in 
lifting  problems  in  the  form  of  a  purely  circulatory  flow.  In  the  study 
of  subsonic  symmetrical  profiles,  this  arbitrariness  In  the  solution 
occurs  when  the  geometry  of  the  wing  is  to  be  determined  from  the  dis¬ 
tribution  of  pressure  exerted  by  the  fluid.  Since,  however,  closure  of 
the  wing  Is  necessary,  an  additional  condition  is  given  which  estab¬ 
lishes  uniqueness  in  much  the  same  manner  that  the  Kutta  condition 
imposes  uniqueness  in  the  lifting  case.  In  the  consideration  of  super¬ 
sonic  quasi— conical  flow,  similar  conditions  to  determine  uniqueness  do 
not  necessarily  arise.  Mathematically,  the  condition  of  uniqueness  is 
determined  from  the  degree  of  dependence  between  simultaneous  linear 
equations. 

One  further  remark  concerning  the  analogy  between  lifting  and  non- 
lifting  problems  appears  to  be  pertinent.  In  two-dimensional  subsonic 
theory  the  integral  relation  between  the  perturbation  velocities  along 
the  wing  surface  may  be  written.  In  the  lifting  case,  as 

w(x,0)  =  ±  rc 

It  Jo  X-Xl 

where  c  is  chord  length.  In  the  symmetrical  case,  the  interrelation 
is  expressed  in  the  form 

u(x,o)  =2.  fC 

Jo 


X— Xi 
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Thus,  aside  from  the  factor  p,  a  complete  duality  exists  in  the  formal 
mathematical  analysis  of  the  two  problems .  Hence,  the  circulatory  motion 
associated  with  a  flat  plate  at  zero  angle  of  attack  (see  sketch)  is 


Lifting  case  Thickness  case 

analogous  to  a  slope  distribution  associated  with  zero  pressure  coeffi¬ 
cient.  However,  in  three-diirensional  supersonic  wing  theory  the  basic 
relations  between  w  and  u  on  the  wing  do  not  have  the  property  of 
duality  (see  equations  (10)  and  (27)).  The  assumption  of  quasi-conical 
flow,  however,  together  with  the  restriction  to  triangular— type  plan 
forms  with  subsonic  leading  edges,  brings  the  study  of  lift  and  thick¬ 
ness  into  more  general  parallelism  and  a  close  similarity  exists 
between  the  final  expressions  in  equations  (25)  and  (32). 

Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  July  10,  1951 
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APPENDIX 

THE  GENERALIZED  PRINCIPAL  PART  AND  FINITE  PART  OF  AN  INTEGRAL 


If  the  singular  integrand  of  a  convergent  improper  integral  is  dif¬ 
ferentiated  formally,  without  due  regard  for  the  singularity,  the  result¬ 
ing  expression  is,  in  general,  improper.  In  applied  theory,  however, 
the  differentiation  is  usually  to  be  performed  upon  the  integral  itself 
and  in  this  case  a  careful  treatment  of  the  entire  expression  leads  to 
a  finite  answer.  The  two  most  common  examples  of  such  problems  arise 
In  the  evaluation  of  the  Cauchy  principal  part  and  Hadamard's  finite 
part.  The  following  development  indicates  the  manner  in  which  these 
cases  are  extended  to  include  multiple  differentiations.  (The  gener¬ 
alized  principal  part  shall  be  concerned  with  integrands  having  singu¬ 
larities  within  the  region  of  integration  and  of  order  n  where  n  is 
a  positive  integer;  the  finite  part,  on  the  other  hand,  involves  inte¬ 
grands  with  singularities  at  an  end  point  of  the  region  of  integration 
and  of  order  n  +  l/2.) 


Consider  first  the  evaluation  of  Cauchy's  principal  part.  In  this 
case  a  single  differentiation  is  used  and  the  expression 


g(x)  = 


_S_ 


A(xi)  In  |x— xi|  cbcij  a<x<b 


■becomes,  for  constants  a  and  b, 

g(x)  =  f*  A<Xl)dXl 
Ja  x“Xl 


(Al) 


(A2) 


Here  the  symbol  on  the  integral  sign  indicates  that  g(x)  is  to  be 
evaluated  by  a  limiting  process  defined  as  follows 


g(x)  = 


lim 

r ~  A(x>>  ^  ♦ 

f*  A(xi)  tax] 

e  ->o 

Lj  a  X-Xl 

'J  x+e  x  Xl 

(A3) 


To  assure  the  convergence  of  this  integral  it  is  sufficient  but  not 
necessary  to  assume  that  A(x)  is  differentiable  at  the  point  xi  =  x 
and  that  elsewhere  within  the  region  of  integration  A(xi)  Is  either 
continuous  or  possesses  integrable  singularities.  The  concept  of  the 
Cauchy  principal  part  is  so  well  known  that  the  symbol  on  the  integral 
is  often  omitted,  as  shall  be  done  here. 
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Turning  next  to  the  case  of  multiple  differentiations,  consider  the 
expressions 


xx— x 


p^°  d  dxx  p  dxx 

*  J  t  a< x< b  <All) 


where  the  limits  of  integration  are  independent  of  x  and  the  symbols 
on  the  two  latter  integrals  indicate  that  the  generalized  principal  part 
is  to  be  calculated.  From  the  definition  of  I  in  the  first  integral 
it  follows  that  * 

I  =  £-  in  tE  =  _  - *=£ - -  (a5) 

x-a  (b— x)  (x-e)  '  } 

The  simple  definition  given  by  equation  (A^)  can  be  generalized  to 
include  integrals  of  the  type 


To  = 


A(zx)  In  jz-Zi]  dzx  =  — 


A(xx)  dxx 


Zi— z 


A(xx)  dxi 
(xx-x)2 


Equation  (A 6)  defines  the  symbol  appearing  in  the  final  member.  It  is 
possible,  however,  to  relate  this  integration  to  the  particular  integrand 
in  the  final  integral  of  equation  (A4)  by  writing  I2  in  the  form 


To  = 


-■*(*>  dIi  +  Mz)  r'jxi 


Zr X 


Xi-z 


Then  if  A(xi)  is  integrable  and  if ,  at  =  x,  its  derivative  exists 
and  is  single  valued,  the  expression  for  I2  becomes 

Ia=  Afa>-  *(x)  toi,A(l) 

JSL  (xi-x)  Ja  (xx-x) 


A(xi)  -  A(x)  ,  A(x)  (b-a) 
- - — ~ —  ——————— 

( xx— x )  (b-x)  (x-a) 


where  the  results  of  equation  (A5)  have  been  used. 
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The  first  integral  in  equation  (A7)  is  now  in  a  form  that  involves 
no  extension  beyond  the  concept  of  Cauchy's  principal  part  and  the 
evaluation  of  I2  may  be  carried  out  with  that  form.  Furthermore,  it 
can  be  shown  from  equation  (A7)  that  if  the  indefinite  integral  of 
A(zi)/(zi— x)2  exists  such  that 


A(x1)dx1 

(xx-x)2 


=  G(xi,x)  +  C 


then  the  value  of  I2  can  be  found  by  following  the  conventional  rules 
for  substitution  of  limits  so  that 


/  =  G(b,x)  -  G(a,x) 

Ja  (xx-x)2 


The  extension  to  higher  ordered  derivatives  is  obvious.  Thus^,  for 
a  and  b  independent  of  x,  one  has 


Sxn+1 


A(xx)  Zn  |x-i xx|  dxx  = 


A(xx) dxx 
(xx-x) 


A(xi)dxx 
x  n+i 


a  (xi-x) 


(A10) 


Equation  (A10)  defines  the  final  member  appearing  in  it.  The  quantity 
In+X  can  also  be  written  in  the  form 


-*-n+i  n* 


A(^i)  ~  B(x,xx) 

/  \n+1 

(xx-x) 


5^dxx  + 


C  B(x,xx)dxx 

„  /  xn+i 

a  (xx-x) 


(All) 


where 


B(x,xx)  =  A(x)  +  (xx-x)  +  .  .  .  +  — - (xx-x) 

1  (n-D! 


r*  gfi  = j_  fb  'f  r^dxi 

a  (xx-x)1+1  ’  il  Z1"X 


7  — “ — r  +  (-!)"  ""  — “ — r  ,i<  n 
*  (b— x)1  (x-a)  J 


b2 
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The  first  n  derivatives  of  A(xx)  are  assumed  to  exist  and  be  single 
valued  at  xx  =  x  while  elsewhere  in  the  range  of  integration  A(x) 

may  possess  integrable  singularities.  The  generalization  of  equation  (A9) 
holds  so  that  if 


♦  c 


then 


(xi-x) 


a  (xx-x) 


(A12) 


(A13) 


It  is  also  possible  to  extend  the  definition  of  equation  (A10)  to 
include  a  functional  dependency  on  x  in  the  numerator  of  the  integrands. 
Thus,  replacing  A(xx)  by  A(x,xx),  equation  (A10)  again  defines  uniquely 
a  principal  part  integral  provided  the  first  n  derivatives  of  A(x,Xi) 
with  respect  to  x  and  x*  exist  at  xi  =  x. 

The  original  concept  of  the  finite  part  was  used  by  Hadamard  in 
connection  with  square  root  singularities.  Consider  the  expressions 


^  ^ a  v  b— Xi 


-f  - 

v.  Jn  Sb 


d  dxi 


0>-*i)3/2 


(Alb) 


From  the  first  integral  in  this  relation  it  follows  that 


The  natural  extension  of  this  idea  is  to  consider 


(A15) 


A(xx)dXi  5  r  a  A(x1)dx1  _  1  ^A(x1)dx1 

^  J-b-Xi  2  i  (b-ijf72 


OD  a  v'b-Xi  a  °b  *  -a 

where  A(xx)  is  continuous  at  xi  =  b  and  is  integrable  elsewhere  in 
the  range  of  integration. 

The  evaluation  of  J2  can  be  related  more  closely  to  the  integral 
of  equation  (Alh)  after  rewriting  J2  in  the  form 
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it  follows  that 


A(xx)dxx 


a 


A(xi)  -  B(b,n) 


'  dxa.  +  J.  -  '  ^  "  (A20) 

a  (b-xx)n4i  Ja  (b-Z!)114^ 


B(h,zi)dzi 


where 


B(bjZi)  =  A(b)  -  A'(h)(b-zi)  +  ...  +  (b-xx)11""1 

(n-l)i 


(f*  d-Xi  _  (—1)^2^  f  —  ^  dzi  _  _  2 _ 1_ 

a  (b-Xi)1+3  1.3. ».(2i— 1)  \dh/  J&  (b-nji  21-1  (b-a)1' 


It  is  furthermore  possible  to  extend  the  definition  of  equation 
(A19)  to  Include  a  functional  dependence  on  b  of  the  integrand  A(xx) 
Replacing  A(xi)  by  A(b,xi),  equation  (A19)  again  defines  uniquely  a 
finite  part  integral  provided  that 

lim  ^(b^xx)  ^ - 

zx  — >b  TTn -  ,  =  0 
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